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One of the necessary covariant conditions for gravitational radiation is the vanishing of the diver- 
gence of the magnetic Weyl tensor H a b, while H a b itself is nonzero. We complete a recent analysis 
by showing that in irrotational dust spacetimes, the condition div H — evolves consistently in the 
exact nonlinear theory. 



Irrotational dust spacetimes, typically considered as models for the late universe or for gravitational collapse, are 
covariantly characterized by the dust four-velocity u a , energy density p , expansion and shear a a b , and by the 
free gravitational field, described by the electric and magnetic parts of the Weyl tensor C a bcd '■ 



E a b = C ac bdU C U d , H a b — hs acc [C' 



cd 



beU 



where h a b = g a b + u a Ub is the spatial projector, g a b is the metric tensor, and e a bc = VabcdU d is the spatial projection 
of the spacetime permutation tensor rj a bcd [0 ■ Gravitational radiation is covariantly described by the nonlocal fields 
E a b , the tidal part of the curvature which generalizes the Newtonian tidal tensor, and H a b , which has no Newtonian 
analogue As such, H a b may be considered as the true gravity wave tensor, since there is no gravitational radiation 
in Newtonian theory. However, as in electromagnetic theory, gravity waves are characterized by H a b and E a b , where 
both are divergence-free but neither is curl-free ||, 

In [0, it was shown that in the generic case, i.e., without imposing any divergence-free conditions, the covariant 
constraint equations evolve consistently with the covariant propagation equations. These equations are: 
Propagation equations 



Constraint equations 



P + e P = o 1 (i) 

<d + \e 2 + \p = -o ab o- ab , (2) 

0~ab + §0CTa6 + E a b = — a^O^ , (3) 

E ab + <dE ab - CUrl H ab + \p° \ib = 3(T c ( a -E b ) C , (4) 

H ab + OH a b + curl E ab = 3<j c{a H b) c . (5) 

D b a ab - |D a 6 = , (6) 

curl o-ab - H ab = , (7) 

D b E ab - iD Q(0 - e abc o- b dH cd , (8) 

B b H ab = -e abc cr b d E cd , (9) 

where S^ a b) — ha c hb d S^ ca ^ — \S c dh cd h a b is the projected, symmetric and trace-free part of S a b , the covariant spatial 

derivative is denned by D a S b "... c = h a p h b q ■ ■ ■ hJV p S q "' ... r , the covariant spatial divergence is D b S a b , and the 
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covariant spatial curl is curl5' a = e a bcD b S c for vectors and curl5 a h = £ c d(aD c Sb) d for tensors. (Further details are 
given in jy, [0.) In the linearized theory of covariant perturbations about an FRW background, the right sides of 
these equations are all zero. 

It was previously claimed that in the exact nonlinear theory, the gravity wave condition 

B b H ab = (10) 

implies H a b = As shown in this claim arises from a sign error and is incorrect, Bianchi type V spacetimes 
providing a counterexample. Here we complete the analysis of [jlj by showing that consistency is maintained if (10) is 
imposed, without H a b zero. 

The fact that consistency is not automatic is illustrated by the case of silent universes, in which H a b = . For 
these solutions, consistent evolution of the condition H a b = imposes a series of nontrivial integrability conditions, 
which are identically satisfied in the linearized case, but not in the nonlinear case @ , B • Thus there is a linearization 
instability in silent universes. By contrast, when (|10|) holds but H a b is not forced to vanish, which includes gravity 
wave solutions, there is no linearization instability following from the evolution of (|To|). An example of consistency 
conditions arising already at the linearized level is given by purely magnetic spacetimes, E a b = , for which Qu a b = 

The proof that ([l0]) evolves consistently is based on a combination of tetrad methods || and the covariant 
methods of JjJ. The only direct effect of ([Io|) on the covariant propagation and constraint equations is an algebraic 
modification of the constraint (|J), which does not change the consistent evolution of the constraints. We have to 
check only consistent evolution of the new condition ( |Io| ) itself. It is more convenient to replace (|l(]) by the equivalent 
condition that follows from (||), 

[<x,E}=0, (11) 



where we are using index-free notation for the covariant commutator. In the linearized case, (|llj) is identically satisfied 
since the left side is second order of smallness, and consistency is automatic. 

In the exact nonlinear case, using only the shear propagation equation (^) and its covariant time derivative, we find 
that 



[a, E] = -[a, a] + §9[ct, E] - a[a, E] 
[&,E] = -ZQ[a > E}+a[a,E}. 

[a,E\ =-[a,a\. (12) 



and 



Adding these equations gives 



Now the right side may be shown to vanish identically without differentiating (|1 1|) , i.e. using only the algebraic 
content of (pi]), as follows. 

From the shear propagation equation (||), ( |ll| ) is equivalent to 

[a,a} = 0. (13) 

We choose an orthonormal tetrad |n]] {eo = u , e^} , with {e M } a shear eigenframe, so that 

ffoa = = d <T 0a , a^y = = <9oo>„ if ^ ^ v , ( 14 ) 

where do denotes the directional derivative along eo = u . Then we have 

[a, &} ab = (<T aa - <Jbb) &ab (nO Slim) . (15) 



At all points where the shear is nondegenerate, i.e., where a aa 7^ Obb when a ^ b, ( p_5[ ) and ( |13[ ) show that a a b is 
diagonal - and thus E a b is also diagonal, by (||). In fact diagonality still holds at points of degeneracy, as follows 
from the tetrad form of the covariant derivative: 

(Tab = do<7ab ~ ^Ob^ac ~ F C o a <7 c b , 

where the Ricci rotation coefficients are T a b c = e a ■ Vf,e c = — T c b a ■ Using ( |l4| ) and (j|), we get 
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a ^ b => <T a6 = (o- oa - (7 fcb ) r b0a = --B ah (no sum) , 



(16) 



so that (T ab is diagonal also where a aa = Obb ( a ^ b). Thus the shear eigenframe simultaneously diagonalizes a a b , 
a a b and E a b . This regains a result given in . 
It also follows from Q, (|f) and © that 

r QOh - o (17) 

holds at all points where the shear is nondegenerate. (Note that (|l7j ) is an identity for a = b .) At points of degeneracy, 
i.e., where an — a 22 > we can use the remaining tetrad freedom of a rotation in the lei, plane to set Tio2 = , 
so that (|l7|) still holds. Specifically, such a rotation through an angle a preserves ( [hf ) and the degeneracy, while 

Ti02 — * Tl02 — d^a ■ 

Thus we can ensure that ( fl^ ) holds throughout spacetime, by specializing the eigenframe where necessary. Then ( |T7j ) 
shows that d a b is also diagonal in this frame, since 

a^b => a ab = (& aa - &bb)T b oa = (no sum) , 

where we have used the fact that do<j a b is diagonal. The covariant (frame-independent) consequence of the simulta- 
neous diagonalizability of a a b and a a b is 

[a, a] =0, 

which shows that the right side of ( p^ ) does indeed vanish identically, consistent with and independent of the derivative 
of ([TT|). Thus the the first covariant time derivative of the condition ( fll|) imposes no consistency conditions. It is clear 
from the above argument that all the subsequent covariant time derivatives of o~ ab are also diagonal in the eigenframe, 
so that these higher derivatives all commute with the shear and amongst themselves. It follows that the second and 
higher covariant time derivatives of the condition ( fill) also vanish without further conditions. 

This establishes that the covariant condition div H = evolves consistently in the exact nonlinear theory. The 
question whether such consistency extends to the further covariant gravity wave condition divi? — is more difficult, 
and under investigation. 

Finally, we note that, by virtue of (|l7|) and the propagation equation (^J), curl H a b is also diagonal in the eigenframe 
that diagonalizes o~ a b and E a b , i.e., there is a shear eigenframe such that a ab , E a b , curl H a b and all their covariant 
time derivatives are diagonal, and therefore commute. 
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